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1 Introduction 



The starting point for most advanced calculations of nonleptonic kaon decays is an 
effective weak lagrangian of the form 0, |^ 

Cl{\AS\ = 1) = V2GfV^,V:,Y.C^O. (1) 

i 

which can be derived with the help of the Wilson operator product expansion from 
elementary quark processes, with additional gluon exchanges. In the framework of 
perturbative QCD the coefficients Cj are to be understood as scale- and normalization 
scheme dependent functions. There exist extensive next-to-leading order (NLO) 
calculations 0, |^ in the context of kaon decays, among others. These calculations 
are based on the possibility of factorization of short- and long-range contributions, 
i.e. into Wilson coefficient functions Ci and mesonic matrix elements of four-quark 
operators Oi, respectively. The latter, however, can presently be obtained only 
by using nonperturbative, i.e. model-dependent, methods. But still it seems not 
possible so far, to fulfil the obvious requirement of scale- and renormalization scheme 
invariance of the resulting amplitudes in a satisfactory way. 

Usually, the results of calculations are displayed with the help of i?-factors in the 
form 



< nnlOAK > 



vac. sat. ) 



where the mesonic matrix elements of four-quark operators are approximated by their 
vacuum saturation values which are, of course, /x- independent. In principle, Bi{fi) 
should be estimated by some higher-order calculations in the long-range regime, 
for instance, in 1 /iVc-expansion in the form 1 + 0{1/Nc), or from the lattice 
approach. The preliminary stage of these calculations is best characterized by the 
problem to explain the well-known A/ = 1/2 rule quantitatively. Of course, the lack 
of such calculations for long-range effects severely restricts the predictive power of 
(|l]), leaving only the possibility of some semi-phenomenological treatment 0, 0, 0, 
with correspondingly large theoretical uncertainties. As a matter of fact, some of 
the combinations [Ci{fi)Bi{fi)] have to be fixed by experimental data, while others 
are restricted by additional theoretical arguments. 

A logical improvement of this approach, pursued in H, P], is to take the lagrangian 
(|lD as defining the structure of the \AS\ = 1 weak interaction and to submit it to 
a global confrontation with data by sandwiching it between appropriate initial and 
final states (for instance, charged and neutral kaons, respectively 27r, Svr, titcj,... 
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states). Of course, there would be hardly any gain in predictive power, if we stuck to 
the vacuum saturation approximation because there may be large unknown factors 
involved, when switching between different initial/final states of the same operator. 
It is evident that everything could be absorbed into the i?j-factors which thereby 
would not only be scheme- and scale-dependent, but would also change with the 
initial/final states (becoming even unknown (complex) functions of dynamic variables 
for final states). On the other hand, these latter changes should be avoided or, at 
least, diminished when using a consistent higher-order calculation of the mesonic 
matrix elements in the framework of chiral theory which today is, certainly, the most 
developed approach to long-range hadronic phenomena in general. These higher- 
order calculations after the inclusion of meson loops have their own scale dependent 
renormalization ambiguities. In the ideal case these meson loops should cancel the 
scale dependence of the Wilson coefficients from the perturbative QCD calculation 
(resulting in 5-factors equal to 1). As a complete match is lacking so far, the best 
one can do is to fix the long-range ambiguities in a definite way at an explicitly 
or implicitly given scale while considering the Wilson coefficients times (unknown) 
i?-factors as phenomenological constants. 

As the above-mentioned perturbative QCD calculations in NLO are believed to 
be reliable down to a scale of 0(1 GeV) (in |^, ^ the charm quark mass rric = 1.3 GeV 
is used as a matching scale), it is very desirable, however, that the effective scale for 
the long-range calculation reaches this value too. Therefore, our strategy includes 
a consistent calculation of the matrix elements up to 0{p^), and furthermore some 
consideration of the effects from vector, axial- vector and scalar resonances which can 



be included by the procedure of reduction [|iy, i.e. a certain recalculation of 
structure coefficients of the effective chiral lagrangians involved. 

In this note we reconsider in the above manner especially the K —>■ 271 chan- 
nels with matrix elements of four-quark operators calculated at 0{p^) in momentum 
expansions within the chiral lagrangian approach. This approach is based on the 
bosonized version of the weak lagrangian (|I|) and chiral effective meson lagrangians 
of higher orders with the structure coefficients fixed by bosonization of the four- 



quark NJL-type interaction. According to Weinberg's power-counting scheme [|T2 



the calculation involves tree-level, one- and two-loop diagrams. The method of su- 
perpropagator regularization [|T^ was used to fix UV divergences arising at the loop 
level. The aim of this note is to give the first numerical results of this calculation 
concerning K ^ 2t[ amplitudes as they lead to new estimates of e' /e (see section 3). 
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In section 2 and appendices A and B we repeat all relevant definitions taken from 
earlier work. 



2 Lagrangians and currents 

In the present paper we use the operators Oi in the representation given in P, 

O2 = ulI^ul dL'j^SL + dil^^UL ulI^sl + 2dLl^dL dil^SL + 2slJ^sl dil^SL , 

O5 = dL-ft^XcSll QR^'^KIr] , 06 = dL-ff,SLi J2 QRl'^lR 

\q=u,d,s j \q=u,d,s 

\q=u,d,s j \q=u,d,! 

where ql^r. = | (1 =F 75)?; K ^^e generators of the SU{Nc) color group; Q is the 
matrix of electric quark charges. The operators O^^q containing right-handed currents 
are generated by gluonic penguin diagrams and the analogous operators g arise 
from electromagnetic penguin diagrams. The operators Oi, 2,3,5,6 and O4 describe 
the transitions with AI = 1/2 and AI = 3/2, respectively, while the operators Oj^s 
contribute to the transition with both AI = 1/2 and AI = 3/2. 

The Wilson coefficients Ci of the effective weak lagrangian (|l|) with four-quark 
operators Oi are connected with the Wilson coefficients Cj corresponding to the basis 
of four-quark operators Qi given in Refs. ll|, §], by the following linear relations: 

Cl = Cl - C2 + C3 - C4 + Cg - Cio , C2 = 7(ci + C2 - Cg - Cio) + C3 + C4 , 

5 

1^ 1/2 



C3 = -C4 = - (^-(Ci + C2) + Cg + Cio^ 



<^5 = C6, C6 = 2(^C5 + ^C6), Cr = ^(^cj + 2csy Cs = ^C8. (2) 

The bosonized version of the effective Lagrangian (|^) can be expressed in the 
form MM: 



= (-6+6+6 



+(6 + 5 6)^ JlMI+ + 4,(4,+ ^J] 
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+^4 
-4^5 



'A,- A,){Jt,+ A,) + 2 JUJl,+ tJl) 



(Jr- A)iJl+ ^4) - iJ] 

2" 



V3 



Jl){Jl+^Jl) 



iJl){V2Jl 



-16^8 



+^^{Jl+^Jl){V2Jl-J\ 



h.c. 



— /8 ) 



^Jl) 



(3) 



Here Gjt' = \/2GFVudV*g, Jl/r^i ^"^^ ^l/r bosonized {V =F A) and {S =F P) 
meson currents corresponding to the quark currents ^7^|(1 =F 7^)A"g and densities 
q\{l =F 7^)A"g, respectively (A" are the generators of the U{?>)f flavor group); 



6 = 

^5,8 



cAi- 



1 



^2,3,4 
1 



2,3,4 



1 + 



2Nr 



^6,7 : 



^6,7 — Cqj 



(4) 



where the color factor 1/Nc originates from the Fierz-transformed contribution to 
the nonleptonic weak effective chiral Lagrangian 

Only the even-intrinsic-parity sector of the chiral lagrangian is required to de- 
scribe nonleptonic kaon decays up to and including O(p^). In the bosonization 
approach this sector is obtained from the modulus of the logarithm of the quark 



determinant of the NJL-type models [15] using the path-integral technique (see [Tl 
and references therein). The meson currents JI/fiij. ^^"^ "^l/r obtained from this 
quark determinant by variation over additional external sources associated with cor- 
responding quark currents and densities P, P]. From the momentum expansion of 
the quark determinant to 0(p^") one can derive the strong lagrangian for mesons 
Ceff of the same order and the corresponding currents J^/RfM and to the order 
0(p^"~^) and 0(p^"~^), respectively. 

For example, from the terms of quark determinant of 0{p^) one gets 



C 



eff - 



-^tT{Ll) + ^tT{xU^ + Ux^] 



^tr(AU 



{pO)a 



F2 

-^mR tr (A"f/) 



(5) 
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where U = exp (^^^(p^ , with (p being the pseudoscalar meson matrix, and = 
D^UU\ D^U = d^U + {Aj:U - UAj^) and Aj^/^ = V,, ± A^ are right /left-handed 
combinations of vector and axial-vector fields. Furthermore, Fq ^ 90 MeV is the 
bare coupling constant of pion decay, x = diag(Xu, X% xi) = — 2mo<gg>-fo~^ is the 

0.47 GeV^ mo is 



meson mass matrix, XZ 



0.0114 GeV^, xl 



0.025 GeV^, xl 



the current quark mass matrix, <qq> is the quark condensate, m 
average constituent quark mass, and R =<qq> / {jtiFq). 
At 0{p^) one gets 



265 MeV is an 



tiL 



+ tr (^-[L^, L,]' + 3{Lir ) + L,ti [{L 



If] 



-U tr {LD tr {xU^ + Ux^) - L, tr (l^xU^ + Ux^ 



+LstT{{x^Uf 
itilx 



{xU^)')+H2tTXX^ 



— m 



UL,ii{xU^ + Ux^) + -L5{L^, {xU^ + t/x^)} 
i?tr {A" [L4f/tr(Lj) + L5(LJ?7) - 2L8f/x^t/ - H^x] } , 



(6) 



where Lj and are structure constants introduced by Gasser and Leutwyler [|T7 



For the sake of brevity, here and in following expressions for the lagrangian and 
currents generated at O(p^), given in appendix A, we restrict ourselves to the terms 
which are necessary to calculate the decay i^' — * 27r at O(p^). We do not show 
explicitly the terms of the effective action at O(p^) generating the scalar current 
Jf^ which is necessary for the full calculation of the tree-level matrix elements at 
O(p^) for the penguin operators, since the corresponding contributions turn out to 
be negligibly small. 

The structure constants Lj, Hi and Qi should be obtained from the modulus of 
the logarithm of quark determinant of the NJL-type model which explicitly contains, 
apart from the pseudoscalar Goldstone bosons, also scalar, vector and axial-vector 
resonances as dynamic degrees of freedom. However, in order to avoid double count- 
ing in calculating pseudoscalar meson amplitudes when taking into account resonance 
degrees of freedom, one has to integrate out (reduce) these resonances in the generat- 
ing functional of the bosonization approach. As a consequence of this procedure, the 
structure coefficients of pseudoscalar low-energy interactions will be quite strongly 
modified. In this way one effectively takes into account resonance-exchange contri- 



butions |jTO|, |TT|. The explicit expressions for the structure constants of non-reduced 
and reduced effective meson lagrangians are given in appendix B. 

In the context of the power counting rules we have to give some comments con- 
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cerning the additional symmetry breaking term ~ tr(mo-D^t/) which appears in 



TU|. This term leads to nonzero meson matrix element of the gluonic penguin 



operator due to the appearance of an additional contribution to the scalar density 
~ tr(D^f/). In the bosonization approach it arises at 0{p^) from that term of the 
divergent part of the quark determinant which also generates the kinetic term of 



div 



167r2 



yti D^{mU + mo) D^{mU + itlq 



X 



(7) 



where D^^W = d^^U^ + {A^W - WA 



m)' y 



A2 

X 



rn^/x with 



X 



<qq> 



—mFQ/{2<qq>). For the phenomenological value of the quark condensate 
-220 MeV one obtains ^ 839 MeV (to be compared with ^ 1020 



1/3, 



MeV in |]). The term ~ tr(mo/^^t/) does not appear explicitly in the effective 
lagrangian (^) because after transformations of double derivatives according to the 
equations of motion 



-D^U D^U^ ■ U - \{Ux^U -x)-\u tr {xU^ - Ux 



2 o 



it is transformed into a combination of terms contributing to the structure coefficients 
L5, Lg and H2 (compare Eq. (|^)): 



^^trix^ D'U + xD'u^) 

X 



167r2 



-xytilLlixU^ + Ux^ 



+ -xytT{ix^Uf + xU^T 
~xy (tr ixU^ -Ux^))' 



xy tr XX 



(8) 



The last term in Eq. (|^) does not contribute in K ^ 27t decays, therefore the corre- 
sponding term (~ L7) in (|^) is dropped. The corresponding numerical contributions 
to the structure coefficients are L5 = 2.9 ■ 10~^, Ls = 1.4 ■ 10~^ and H2 = —2.9 • 10"^. 
The complete expressions (|T5|) in appendix B contain further contributions to L5, 
Lg and H2, which originate from the divergent as well as the finite parts of the 
quark determinant at O(p^), and result in the numerical values L5 = 0.98 ■ 10"'^, 
Ls = 0.36 ■ 10-3 and H2 = 1.01 ■ lO'^ (see table 6 in appendix B). 
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3 Amplitudes and phenomenological results 



Using isospin relations, the K ^ 2tx decay amplitudes can be parameterized as 



V3 . 



■An 



The isotopic amplitudes ^2,0 determine the K ^ 271 transitions into states with 
isospin J = 2, 0, respectively: 



A, 



02 e 



id2 



Ao = ao e 



where ^2,0 ^'I'e the phases of Trvr-scattering. It is well known that direct CP violation 
results in an additional (small) relative phase between a2 and oq. Let us next in- 
troduce the contributions of the four-quark operators Oi to the isotopic amplitudes 



W'P by the relations 



Ai 



(9) 



1=1 



where ^2 = V2J^o = ^GpFoiml - ml). 

At 0{p^) we obtain for the nonzero tree- level amplitudes A^f* the following ex- 
pressions: 



4^^ 
aS^ 



Af 



__4{2,3) . 



> 



m 



K 



rat 



m\ — ml 



^0 



,2 I ,,2 



QL^ixi + xi + xl) 

+ {L, - ALs)ixl + 3x' + 2x1) + '2L,ml 
QL^ixl + xl + xl) 

+ {L, - 4Ls){xl + 3Xd + 2x1) + 2L,ml 



(10) 



The Ls and H2 contributions in the penguin operators O^^s also have a tadpole 
contribution from K — > (vacuum) included through strong rescattering, K ttttK 
with K (vacuum). At O(p^), in case of the penguin operator O5, the Lg and 
H2 contributions in the direct matrix element from K ^ 27t vertices, are fully 
cancelled by the tadpole diagrams Q. This is due to the possibility to absorb the 
^We thank W.A. Bardeen and A.J. Buras for drawing our attention to this point. 
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tadpole contribution into a redefinition of the — > 27r vertex if all particles are on 
mass shell. In higher orders, however, such a cancellation is generally not expected. 
Therefore, one of the main aims of the present paper is to perform all calculations 
up to 0{p^) to clarify this point quantitatively. 

Some interesting observations on the difference of the momentum behavior of 
penguin and nonpenguin operators can be drawn from power-counting arguments. 
According to Eq. (J^) the leading contributions to the vector and scalar currents are of 
0{p^) and 0(p°), respectively. Since in our approach the nonpenguin operators are 
constructed out of the products of (V^ — A)-currents J£^, while the penguin operators 
are products of (5* — P)-currents J£, the lowest-order contributions of nonpenguin 
and penguin operators are of O(p^) and 0(p°), respectively. However, due to the 
well-known cancellation of the contribution of gluonic penguin operator at the 
lowest order |^, the leading gluonic penguin as well as nonpenguin contributions 
start from O(p^) [|. Consequently, in order to derive the {V — A)-currents which 
contribute to the nonpenguin transition operators at leading order, it is sufficient 
to use the terms of the quark determinant to 0{p^) only. At the same time the 
terms of the quark determinant to 0{p'^) have to be kept for calculating the penguin 
contribution at 0(j9^) since it arises from the combination of {S — P)-currents from 
Eqs. (H) and (|^), which are of OijP) and O(p^), respectively. In this subtle way 
a difference in momentum behaviour is revealed between matrix elements for these 
two types of weak transition operators; it manifests itself more drastically in higher- 
order lagrangians and currents. This fact makes penguins especially sensitive to 
higher order effects. 

Using the truncated lagrangian (|^) only, the gluonic penguin operator matrix 
element at O(p^) is 



However, taking into account the additional contribution from other parts of the 
quark determinant at O(p^), we get a suppression of this matrix element which after 
substitution of the full expression for the non-reduced structure constant L5 from 
Eq. (|15D becomes 



^There is no cancellation of the contribution of the electromagnetic penguin operator Cg at the 
lowest order and the first terms in the expressions ( |lO| ) for A^^\ correspond to the contributions at 



AT =4:R^ -20.0 . 



(11) 
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The modification of the structure constant after reduction of resonances {L^^'^ = 
1.64 ■ 10~^, see appendix B) leads to an increase in absolute value of the gluonic 
penguin matrix element to Aq ~ —11.2, which, however, is still about factor 2 
smaller than (0). 

Table 1 presents the modification of the amplitudes A^p when including suc- 
cessively the higher order corrections O(p^) and 0{p^) and the reduction of meson 
resonances. Our calculations involve Born and one- and two-loop meson diagrams 
and take into account isotopic symmetry breaking (vr^ — 77 — 77' mixing). In our 
approach, the UV divergences resulting from meson loops at O(p^) and 0{p^) were 
separated using the superpropagator regularization method ||13| which is particularly 
well suited to the treatment of loops in nonlinear chiral theories. The result is equiv- 
alent to the dimensional regularization technique, the difference being that the scale 
parameter fi is no longer arbitrary but fixed by the inherent scale of the chiral theory 
/i = AttFq ^ 1 GeV, and the UV divergences have to be replaced by a finite term 
using the substitution 

{C - 1/e) ^ CsP = -I + + pn , 

where C = 0.577 is Euler's constant, e = {A — D)/2 and (3 is an arbitrary constant in- 
troduced by the Sommerfeld- Watson integral representation of the superpropagator. 
The splitting of the decay constants and is used at 0{p^) to fix Csp ~ 3.0. 

The strong interaction phases 62^ arise first at 0{p'^), but for the quantitative 
description of the phases it is necessary to go beyond 0{p^). At O(p^), for the vrvr- 
scattering phase shifts and their difference A = 5o — ^2, we have obtained the values 
of 60 ~ 22°, 62 ~ —13°, A ^ 35° which are in agreement with At O(p^), we have 
obtained 6q ~ 35°, 62 ~ —9°, A ^ 44°, in better agreement with the experimental 
value A^^P = (48 ± 4)° [|2 . 



In our approach the parameters in Eq. (|^) are treated as phenomenological 
(//-independent) parameters to be fixed from the experimental data. They can be 
related to the //-dependent QCD predicted ^i(/t) with the help of some /i-dependent 
i?j-factor defined as 

= 6(/x)i?.(/i). 

Table 2 shows the QCD predictions for the coefficients ^i(/t) = ^i^\lJ') + T^i^\lJ') 
which correspond to the Wilson coefficients 

V V* 
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from the table XVIII of Ref. |I| calculated numerically from perturbative QCD at 
/i = 1 GeV for rrit = 170 GeV in leading (LO) and next-to-leading order in different 
renormalization schemes (NDR and HV). and ^f^^ were obtained from Zi and yi, 
respectively, using the Eqs. and (^). 

As we cannot calculate the factors Bi{fi) theoretically, they can only be fixed 
from data in the spirit of the semi-phenomenological approach JI], ^, 0]. Table 1 
shows that the amplitudes of K 2tt decays are dominated by the contribution of 
the operators Oi with i = 1,2, 3, 4, 5, 8. Moreover, in case of the operators Ci,2,3, the 
first term in the combination (— ^2+^3) dominates in the effective weak meson 
lagrangian (|^). Thus, the isotopic amplitudes can be given in the approximation of 
the dominating contributions of four-quark operators as 



,(1) , A^^y)( ,,\-R ( ,,\ 



+d^'^V)i?5(/x) 4'^ + d^'^V)i?8(/x) 4'^ , (12) 



and 



At least two factors Bi and -B4 can be estimated from the experimental values A'^^ ~ 
10.9 and ^2^*' ~ 0.347 while the other two (penguin) factors B^ and B^ should be 
fixed from other data or restricted by theoretical arguments. 

The parameter e of direct CP- violation in — ^ 27r decays can be expressed by 
the formulae 

uj Imao. r^^i(n/2+5.~5.) . Reaa ^ llmoa 



£ = — - 



\f2 Re ao ' Re oq ' Im oq ' 

and the ratio e je can be estimated as 

{y) 

- = lm\t{Po-P2), Pj = ^^ (13) 

with ImAj = lmV*gVtd = \Vts\\Vtd\sin5 in the standard parameterization of the CKM 
matrix. In our estimates for e' /e we will use the restriction ^ 

0.86 ■ 10"^ < Im At < 1.71 ■ 10"^ (14) 

obtained from the phenomeno logical analysis of indirect CP violation in i^T —> 27r 
decay and B^ — 5^ mixing. 
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Table 3 gives the estimates of e je from a semi-phenomenological approach ob- 
tained after fixing the correction factors B\ and i?4 for isotopic amphtudes in the 
representation ([T2|) by experimental (CP-conserving) data on Re ^0,2, and setting 

= = \. We have used the matrix elements of the operators Oi displayed 
in table 1, and the theoretical values .^j(yu) from table 2. Taking into account the 
dependence of the results on the renormalization scheme, we have obtained without 
reducing resonances (see table la) the following upper and lower bounds for e je^ 
corresponding to the interval (|14D for ImA^: 

-4.7 ■ 10"^ < e je < -0.4 ■ 10"^ . 

The peculiarity of these results lies in the observation that all estimates of e' /e 
lead to negative values. This is related to the fact that in the case corresponding 
to table la the contribution of gluonic penguins to AI = 1/2 transitions appears 
to be suppressed, leading, after the interplay between gluonic and electromagnetic 
penguins, to the relation Pq < P2 for the two competing terms in ([T3|). Generally 
speaking, AI = 1/2 transitions loose importance compared to AI = 3/2 when 
estimating e' /e. The situation changes after the reduction of resonances, due to a 
relative enhancement of the matrix element for the operator O5 (see table lb). The 
bounds obtained in this case are 

-3.0 • 10"^ < e'/^ < 3.6 ■ 10-\ 

that substantially agrees with the bounds given by [|I|, 0]. 

So, our calculations have shown that especially the penguin matrix elements 
are most sensitive to various refinements: higher-order derivative terms in chiral 
lagrangians, the reduction of meson resonances, 7r° — 77 — 77' mixing, meson loop 
corrections. It should be added that the modification of penguin matrix elements, 
discussed in this note, is much more important for gluonic than for electromagnetic 
penguin transitions. This is obvious from the observation that the latter at lowest 
order contain terms of 0{p^) which are left unchanged when taking into account the 
additional terms derived from the quark determinant at O(p^). 

Finally, we give some results concerning the dependence of the above semi- 
phenomenological estimates for e' /e on the choice of the penguin correction factors 
i?5 and -Bg- In table 4 we present the i?5-dependence of e'/e estimated with reducing 
meson resonances for the central value ImAj = 1.29 ■ 10""^ and Bg = 1. The cor- 
responding Pg- dependence of e' /e {B^ = 1) is given in table 5. The sensitivity to 
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a small (~ 5%) uncertainty from the phenomenological fixation of <qq>^^^ is also 
shown in figures 1 and 2. It turns out to be problematic to obtain the value of e' /e 



as large as (20 — 40) ■ 10~^ and consistent with the NA31 result |23| even varying the 
input parameters <qq>^^^, and in a wide region. Taking B^ = B^ = 1^ a. value 
of I <qq>^/^ I > 300 MeV would be required in order to reach e /e > 20 ■ 10~^, which 
is far above its phenomenological limits derived within the usual chiral framework. 

Since our results are very sensitive to the relative contribution of the gluonic 
penguin operator, the question of its phenomenological separation in K ^ 2ti de- 
cays becomes critical, in the context of the A/ =1/2 rule as well as for the very 
important problem of direct CP- violation, where at least one experimental result 
would lead to some revision of the present picture. C7P-conserving K ^ 271 data 
alone are clearly not sufficient for such a separation. It could be accomplished, on 
the other hand, when taking into account Dalitz-plot data for i^' ^ Stt as well as dif- 
ferential distributions for radiative decays K 2717, K vr27, which are described 
by the same lagrangian (|1|). As emphasized above, the reason for this possibility 
is found in the difference in momentum power counting behaviour between penguin 
and non-penguin matrix elements, which appears in higher orders of chiral theory, 
when calculating various parameters of differential distributions, for instance, slope 
parameters of the Dalitz-plot for i^' ^ Svr. A substantial improvement in the accu- 
racy of such experimental data (mostly being of older dates) would be very helpful 
for such a phenomenological improvement of the theoretical situation for e' /e (see 
for a discussion of this point and [0, ^ for some recent measurements). Of course, 
for all these model developments the new experiments at CERN, FNAL and planned 
at Frascati have to be considered as crucial if the accuracy level of 10~^ for e' /e is 
obtained. 

The authors gratefully acknowledge fruitful and helpful discussions with W.A. 
Bardeen and A.J. Buras. 
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Appendix A 



rip') 



At O(p^) one gets 

+Qn[x(D^D,U^L^'L'' + R'R^UD^D.U^) 
+x\D^,D,URm'' + UL^D^D.U) 

+Qi^X^L^xR^ + Q^,(x^xR^.R'' + XX^^M^" 



+Qi9[{xU^f + (x^t/)'] + Q2o(f/^xx^x + Ux^x^ 

where Qi are structure constants introduced in whereas R^ = U'^D^U. The 
corresponding terms of {V =F A) and (S* =F P) bosonized meson currents are given by 



Ljj. 



i-tr<| A" 



UWUU^ + x^U)D,D,U^ + D,D,U\Ux^ + xU^)] 



+ L-UUU^ + x^U)D^D,U^ + D^D,U\Ux^ + xf/^) 



+ /^'^(([/x^ + xU^)DuD,U + D,D^f/(f/tx + X^U))U\ 
+2Q,,{Ux^L^xU^ + xU^L^Ux^) + 2QMUx\U\ L^} + {xx^ ^^4) 
+2QM{Ux^)\L,} + {{xU^)\L^}) 
-AQ,siUx^L,Ux^ + xU^L,xU^) 



and 



J 



(p4)a 



mi?tr{A'^[Qi2 L''U{D^,D^}U^ UR"" + Qi^iL'^L'^D^DM + L'^L'^f/ ■ P'^P'^) 

+Qii{UD''WU^D^D^U + D^D^UD'^WU^ ■ U) 

+gi5 + ^?i6 (X^J + ^'x) + QniUx^URl + Lpx^U) 

+2gi8 L^^Ux^L^u + Qi.iUx^fu + g2o(xf/^x + xx^f/ + f/x^x)] } • 



■^The rather lengthy full expression for bosonized effective lagrangian at 0{p^) was presented in 
refs.p^,^^. 
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Appendix B 



Without reduction of resonance degrees of freedom the structure constants Lj = 
iV,/(167r2) ■ li, H2 = NJilGn"^) ■ hi and Qi = NJ {327r'^m'^) ■ are fixed from the 
bosonization of an NJL-type model as 



1 2 ^ 24 6 



xy — X . 



I, 



1 2 1 , 2 1 levr^xF^ 
-xy — X y , /lo = —xy — 2a; y -\ 1 — — ^ , 

2 ^ ^ 24 ' ^ ^12 Nr Am^' 



(15) 



and 



qi2 
qi5 



1 

60 
2 

3 

1 



«?13 

x(l — x) 



3 V20 



a; + c 



gi4 



gi6 



" 6' 

1 42 a: 



120 



+ -x^ + -(1 - 4x) - 2 X - 



6 



c, 



Qn = h -(1 - 4x) 

120 6^ ' 



4 



+ g )c, gi8 = 772; + ( - - X )c 



1 



gi9 

^20 



240 



X 



-x'^ + x(l + 2xy)c . 



240 



+ x^ + 2(1 - 2y)x^ - x(l + 2xy)c , 



where x = —mF^/{2<qq>), y = Att'^Fq /{Ncm'^) = 1.5 and c = 1 — l/(6y). After 
reduction of the resonances, the structure coefficients get the form 



r 



red 



2 ^ 12 



Z\ + 2{Z\-l){U,{Z\-l)-Z'^ 



jred 
'3 



Zi + 3(Zl-l)Qy(Zl-l)-Zl)" 



and 



„red jred n „red ^ ryQ 



red „reii 

yi6 ~ yi7 



qlf = -2qlt = ^Z\ 



4-6(l + 4(l-Z2))(l-y) +4(1 + 16(1-Z 
3y-2Zl(5-12{l-Z^ 



72 I r- ion ^2^(1 



red n^^<i 

Vig ~ 2^18 



192 



Zl(3y - 2) 



where = 47r^FQ^/ {Z\Njff) = 2.4, and = 0.62 is the tc — Ai mixing factor. 

The current experimental status of the effective chiral lagrangian at O(p^) has 
been discussed in some detail in [03. In table 6. we present the predictions of the NJL 
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model for the structure coefficients Li and H2 (without and with reduction of meson 
resonances) in comparison with their phenomenological vahics. The predictions of the 
NJL model after reduction of meson resonances turn out to be in a good agreement 
with phenomenology. 
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Table 1 . Isotopic amplitudes of the K 27: decays 
with successive inclusion of higher-order corrections. 



a) Without reduction of meson resonances: 







Oi 


O2 




Oi 


Or, 




<\0- 




0(p^) 




-1.000 


1.000 


1.000 


0.000 


-6.833 


0.000 


0.016 


1.015 






0.000 


0.000 


0.000 


1.000 


0.019 


0.000 


-0.016 


0.454 






-1.184 


1.173 


1.090 


0.021 


-7.836 


0.003 


0.016 


1.206 






-0.482 


0.482 


0.482 


0.000 


-3.229 


0.000 


0.008 


0.534 






-0.005 


0.016 


0.035 


0.867 


-0.077 


-0.003 


-0.015 


0.424 






0.000 


0.000 


0.000 


-0.213 


-0.004 


0.000 


0.003 


-0.064 




Re^*^ 


-1.014 


1.003 


0.888 


0.022 


-6.413 


0.003 


0.012 


1.042 






-0.707 


0.709 


0.682 


0.000 


-4.582 


-0.001 


0.011 


0.740 




Re4'^ 


-0.004 


0.016 


0.035 


0.812 


-0.078 


-0.003 


-0.015 


0.396 




luiA^2^ 


0.000 


0.001 


0.001 


-0.114 


0.001 


0.000 


0.002 


-0.054 



b) After reduction of resonances: 









O2 


O3 


C4 


0, 


0, 








ReA^^ 


-1.000 


1.000 


1.000 


0.000 


-11.248 


0.000 


0.016 


1.325 




Re^^^^ 


0.000 


0.000 


0.000 


1.000 


0.017 


0.000 


-0.016 


0.574 




ReAS^ 


-1.195 


1.187 


1.101 


0.021 


-13.524 


0.002 


0.016 


1.554 


0{p') 




-0.482 


0.482 


0.482 


0.000 


-5.365 


0.000 


0.008 


0.617 




Re4'^ 


-0.006 


0.014 


0.035 


0.879 


-0.124 


-0.002 


-0.015 


0.502 






0.000 


0.000 


0.000 


-0.213 


-0.004 


0.000 


0.003 


-0.094 




ReA^o^ 


-0.974 


0.967 


0.798 


0.024 


-10.447 


0.002 


0.010 


1.328 






-0.643 


0.646 


0.619 


0.000 


-7.140 


-0.001 


0.010 


0.864 




Re4'^ 


-0.005 


0.014 


0.036 


0.779 


-0.138 


-0.002 


-0.015 


0.479 






0.000 


0.001 


0.001 


-0.126 


-0.007 


0.000 


0.002 


-0.066 
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Table 2. QCD predictions for the parameters 
= + calculated with Wilson coefficients 

Cj(/i) = Zi{y) + ryi{ij) at = 1 GeV for = 170 GeV 





= 215 MeV 


Ag^ = 325 MeV 


= 435 MeV 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


SI 


-1.286 


-1.061 


-1.165 


-1.443 


-1.159 


-1.325 


-1.624 


-1.270 


-1.562 




0.187 


0.195 


0.198 


0.172 


0.176 


0.182 


0.157 


0.150 


0.165 


S3 


0.129 


0.143 


0.137 


0.122 


0.137 


0.130 


0.115 


0.131 


0.121 


S4 


0.645 


0.714 


0.687 


0.609 


0.684 


0.650 


0.573 


0.654 


0.599 


A^) 

S5 


-0.008 


-0.020 


-0.008 


-0.012 


-0.032 


-0.013 


-0.016 


-0.056 


-0.023 


S6 


0.000 


-0.003 


0.000 


-0.001 


-0.007 


-0.001 


-0.002 


-0.021 


-0.007 




0.002 


0.003 


-0.001 


0.004 


0.008 


0.001 


0.006 


0.015 


0.032 




0.000 


0.002 


0.001 


0.001 


0.004 


0.002 


0.001 


0.009 


0.067 


SI 


0.044 


0.038 


0.048 


0.054 


0.048 


0.053 


0.065 


0.060 


0.069 


S2 


-0.028 


-0.029 


-0.030 


-0.029 


-0.033 


-0.030 


-0.030 


-0.033 


-0.030 


Ay) 

S3 


-0.002 


-0.002 


0.001 


-0.002 


-0.002 


-0.002 


-0.002 


-0.002 


-0.002 


S4 


-0.009 


-0.010 


0.004 


-0.008 


-0.009 


-0.009 


-0.008 


-0.009 


-0.008 


S5 


-0.081 


-0.076 


-0.067 


-0.109 


-0.111 


-0.092 


-0.143 


-0.173 


-0.132 


Ay) 

S6 


-0.033 


-0.042 


-0.021 


-0.049 


-0.076 


-0.033 


-0.071 


-0.139 


-0.051 




0.033 


0.004 


0.006 


0.044 


0.013 


0.016 


0.057 


0.027 


0.032 




0.031 


0.028 


0.031 


0.043 


0.041 


0.045 


0.058 


0.061 


0.067 
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Table 3. Predictions for the parameters of ^ 27r decays 
in the semi-phenomenological approach (^5 = Bg = 1). 
The ratio e' /e is given in units 10~^. 



a) Without reduction of meson resonances: 





= 215 MeV 

MS 


= 325 MeV 

MS 


= 435 MeV 

MS 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Ba 
Po 
P2 

(e / S^rnin 
(£ / E^max 


5.52 
0.65 
0.95 
1.71 
-1.3 
-0.7 


6.28 
0.59 
0.45 
1.46 
-1.7 
-0.9 


5.90 

0.62 

-0.09 
2.64 
-4.7 
-2.4 


5.07 

0.69 
2.09 
2.76 
-1.2 
-0.6 


5.91 

0.61 
1.63 
2.65 
-1.7 
-0.9 


5.39 

0.65 
1.12 
2.68 
-2.7 
-1.3 


4.63 

0.74 
3.59 
4.13 
-0.9 
-0.5 


5.52 

0.64 
4.21 
4.65 
-0.8 
-0.4 


4.73 
0.70 
2.96 
4.28 
-2.3 
-1.1 



b) After reduction of resonances: 





a(^ = 215 MeV 


= 325 MeV 


A% = 435 MeV 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


Bi 


5.85 


6.62 


6.25 


5.36 


6.20 


5.69 


4.88 


5.72 


4.97 


Ba 


0.68 


0.61 


0.64 


0.72 


0.64 


0.67 


0.76 


0.66 


0.73 


Po 


3.21 


2.56 


1.79 


5.12 


4.69 


3.72 


7.59 


8.94 


6.67 


P2 


2.77 


2.43 


3.54 


4.21 


4.07 


4.00 


6.06 


6.85 


6.05 


(s / £)min 


0.4 


0.1 


-3.0 


0.8 


0.5 


-0.5 


1.3 


1.8 


0.5 


(e / S)rnax 


0.8 


0.2 


-1.5 


1.6 


1.1 


-0.2 


2.6 


3.6 


1.1 
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Table 4. Dependence of the predictions for the ratio e je 
(in units 10~^) on the factor ^5 {B^ = 1). 





= 215 MeV 


aJL = 325 MeV 

MS 


^WS = MeV 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


1.0 

1.1 

1.2 
1.3 
1.4 
1.5 
1.6 
1.7 
1.8 
1.9 
2.0 


0.56 
1.10 
1.64 
2.17 
2.71 
3.25 
3.78 
4.32 
4.86 
5.39 
5.93 


0.16 
0.67 
1.17 
1.67 
2.17 
2.67 
3.17 
3.68 
4.18 
4.68 
5.18 


-2.26 
-1.81 
-1.37 
-0.92 
-0.48 
-0.03 
0.42 
0.86 
1.31 
1.75 
2.20 


1.17 
1.89 
2.60 
3.32 
4.04 
4.76 
5.47 
6.19 
6.91 
7.63 
8.34 


0.80 
1.52 
2.25 
2.97 
3.70 
4.42 
5.15 
5.87 
6.60 
7.32 
8.05 


-0.36 
0.25 
0.86 
1.47 
2.08 
2.69 
3.30 
3.92 
4.53 
5.14 
5.75 


I. 97 
2.92 
3.86 
4.81 
5.75 
6.69 
7.64 
8.58 
9.52 
10.47 

II. 41 


2.70 
3.82 
4.94 
6.06 
7.18 
8.30 
9.42 
10.53 
11.65 
12.76 
13.88 


0.80 
1.70 
2.59 
3.48 
4.38 
5.27 
6.16 
7.06 
7.95 
8.84 
9.73 



Table 5. Dependence of the predictions for the ratio e /e 

(in units 10^^) on the factor {Br^ = 1). 





= 215 MeV 

MS 


^Ws = MeV 


A^ = 435 MeV 

MS 




LO 


NDR 


HV 


LO 


NDR 


HV 


LO 


NDR 


HV 


1.0 
1.1 
1.2 
1.3 
1.4 
1.5 
1.6 
1.7 
1.8 
1.9 
2.0 


0.56 
0.30 
0.03 
-0.23 
-0.50 
-0.77 
-1.03 
-1.30 
-1.56 
-1.83 
-2.09 


0.16 
-0.07 
-0.31 
-0.55 
-0.78 
-1.02 
-1.25 
-1.49 
-1.73 
-1.96 
-2.20 


-2.26 
-2.52 
-2.78 
-3.05 
-3.31 
-3.57 
-3.84 
-4.10 
-4.36 
-4.62 
-4.89 


1.17 
0.80 
0.43 
0.07 
-0.30 
-0.67 
-1.04 
-1.40 
-1.77 
-2.14 
-2.51 


0.79 
0.45 
0.10 
-0.25 
-0.60 
-0.94 
-1.29 
-1.63 
-1.98 
-2.32 
-2.67 


-0.36 
-0.75 
-1.13 
-1.51 
-1.89 
-2.28 
-2.66 
-3.04 
-3.42 
-3.80 
-4.18 


1.97 
1.48 
0.98 
0.48 
-0.02 
-0.52 
-1.02 
-1.52 
-2.01 
-2.51 
-3.01 


2.70 
2.19 
1.67 
1.15 
0.64 
0.12 
-0.39 
-0.90 
-1.42 
-1.92 
-2.43 


0.80 
0.32 
-0.15 
-0.61 
-1.07 
-1.52 
-1.96 
-2.40 
-2.82 
-3.25 
-3.66 
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Table 6. Theoretical and phenomenological values of the structure 
coefficients Lj and H2 ^ (in units 10^^). 



Structure 


Without reduction 


After reduction 


Phenomenology 25] 


coefficients 


of resonances 


of resonances 




Li 


0.79 


0.85 


0.4 ± 0.3 


L2 


1.58 


1.70 


1.4 ±0.3 


L3 


-3.17 


-4.30 


-3.5 ± 1.1 


L4 








-0.3 ±0.5 


L5 


0.98 


1.64 


1.4 ±0.5 


Ls 


0.36 


1.12 


0.9 ±0.3 


H2 


1.01 


-0.67 





^The parameter H2 has not been determined phenomenologically until now. 
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